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Dielectric responces of the one-dimentional electron system is investigated numerically. 
We treat an interacting one-dimentional spinless fermion model with disorder by using the 
Density Matrix Renormalization Group(DMRG) method which is extended for nonuniform 
systems. We apply an electric field E to the system and calculate dielectric responces. Dielec- 
tric responce of the Mott insulator and the Anderson insulator are calculated respectively. 
Steplike behaviors in P ~ E curve are obtained which corresponds to breakdown of the in- 
sulating behavior. For the Mott insulator, the steps originate from generation of kink-pairs. 
For the Anderson insulator on the other hand, the origin of the steps is a crossing of the 
localized one particle energy levels. We also treat random systems with interaction. From one 
parameter scaling analysis of the susceptibility the metal-insulator transition in attrac- 
tively interacting region is confirmed and a phase diagram of the random spinless fermion 
model is obtained. 

KEYWORDS: Anderson localization, polarization, dielectric response, density matrix renormal- 
ization group 



1. Introduction 

Effects of randomness and interaction in electronic systems are two major problems in con- 
densed matter physics. These two problems have been studied intensively for several decades 
and a lot of fundamental results are accumulated. The presence of disorder entails a local- 
ization of electronic states due to a quantum mechanical interference of the Bloch states and 
this phenomenon has been well studied as the Anderson localization. According to the scaling 
theory, all states localize in one- or two-dimensional systems without interaction, no matter 
how weak the randomness is. Interaction between electrons also changes states of electrons 
drastically. Strong Coulomb interaction can lead a metallic system to an insulator when the 
filling factor is rational and attractive interaction may cause superconductivity. Although we 
have well established understanding of the localization of non-interacting electrons, study on 
the correlated electrons with randomness is still unsatisfactory. It is mainly due to lack of 
reliable numerical techniques for the correlated electrons with randomness. 

Today, for the one-dimensional correlated systems without randomness, a lot of numer- 
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ical results are collected by the Density Matrix Renormalization Group (DMRG)^^ method 
and consistent understanding with analytical predictions has been achieved. "^^ Recently, the 
DMRG was extended to random systems and it has been applied to random spin chains.^) 
It brings us new possibilities to investigate correlation effects on the random system. In this 
paper, we concentrate on dielectric responses of one-dimentional electronic systems with both 
randomness and interaction. Applying a finite electric field E, we calculate a polarization P 
as a derivative of ground state energy by the field strength E. Linear susceptibility x is also 
obtained numerically as a derivative of the polarization by the electric field. The polarization 
P reflects deviation of mass center of electrons by the field. The dielectric response is a fun- 
damental quantity to characterize whether the system is a metal or an insulator. It is suitably 
defined in open boundary conditions, while Resta and Sorella defined polarization in peri- 
odic system recently.^-* Hence the DMRG is naturally applied to investigating the dielectric 
response. Using the conventional DMRG, Aebischer et al. investigated the dielectric responce 
of Hubbard model with next nearest neighbor hopping by and confirmed the metal-insulator 
transition of the system.^) We apply the DMRG for nonuniform system by Hida"^-* to random 
fermionic systems. 

We treat an interacting one-dimensional spinless fermion model in random potential. The 
Hamiltonian of the L site system is given as 

L-l L-l L 

H = -t ^(ctcj+i + h.c.) + V^Y^ riirii+i + ^ eiUi, (1) 

1=1 1=1 i=l 

where rij = c|ci and is a random potential, which distributes over an interval [—W/2, W/2] 
uniformly. We set t = 1 and consider the half-filling case and impose an open boundary 
condition. In the absence of disorder, the system is metallic for —2 < V < 2. For half-filling 
case, &t V = 2 the system undergoes a metal-insulator transition and the system has a finite 
charge gap for V > 2. In an attractive interaction region sit V = —2, the system becomes 
unstable due to phase separation. When the randomness is present, the system is always 
insulating due to the Anderson localization without interaction. Then the interplay between 
the randomness and interaction can be interesting^' and possible metallic phase is expected 
for some range of negative V. Chui and Bray, and Apel treated the effect of interplay between 
randomness and interaction on Tomonaga model. ^' The former authors determined critical 
value of the interaction constant of localize-delocalize transition by analysis of the density 
autocorrelation function. The latter considered the dynamic conductivity. The spin-dependent 
backward scattering was also treated by Chui and Bray,^^) Apel and Rice,^^) and Suzumura 
and Fukuyama.^"^' Gianiarchi and Schulz took into account the renormalization of interaction 
by the disorder and obtained a phase diagram. Runge and Zymanyi, and Bouzerar and 
Poilblanc estimated the size of delocalized region in random interacting spinless fermion model 
by the exact diagonalization.^^'^^^ Schmitteckert et al. used the DMRG and estimated the size 
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of the delocalizcd region from analysis of the phase sensitivity. 

In sec. II, we discuss on the dielectric responses of systems. Sec. Ill is for results and 
discussion. Sec. IV is a summary. 

2. Dielectric response 



We focus on dielectric responses of the one dimensional system. In order to observe the 
dielectric response directly, we apply the electric field E to the system. As a second quantized 
form of the potential, —Ex, the coupling term 

HE = -E^(i-^)en,, (2) 
i=l ^ 

is added to the Hamiltonian, where e is the charge of the electron and we set e = 1 in this 
paper. Then the full Hamiltonian of the system is given by Ht = H + He- As a function of 
E, the polarization P of the system is defined as 

-4^ = -km. 



i=l ^ ' 

where Eq is the ground state energy and (nj) e represents the ground state expectation value 
of rii.^^ Here we used the Feynman's theorem to derive the expression. 

In the Mott insulator, for a finite value of the electric field whose potential energy is 
comparable with the Mott gap, we expect a collapse of the charge gap due to the interaction. 
Then one can expect a gap closing which is a collapse of the Mott gap at Eg ~ EL where Eg is 
the Mott gap. In the Anderson insulator, reconstruction of the charge occurs by transferring 
electrons above the tunneling barrier. The critical field strength Ec is estimated as Ec ~ E^/L 
where E^ is an energy difference between the highest occupied and the lowest empty one 
particle states. E^ is of the order of 1/L on the average. In each case, we can obtain information 
on the charge degree of freedom above its ground state. 

As for a linear response regime, E ^ 0, we calculate zero-field dielectric susceptibility as 



_ dP 

^~ dE 



1 d'^Eo 



E=o L dE^ 



(4) 

E=0 



Prom the susceptibility %, we directly obtain information whether the ground state is 
metallic or an insulator. In the thermodynamic limit, % is diverging if the system is metallic, 
but converging to a finite value if it is an insulator. Indeed % ~ is expected by the per- 
turbation calculation for the pure non-interacting system(W^ = 0, F = 0). It is also confirmed 
numerically for the pure interacting systems(T^ = 0, — 2 < F < 0) by the DMRG(See later). 

In order to calculate the charge distribution and the ground state energy, we use the exact 
diagonalization for V = 0, and the DMRG for the finite V. In the application of the DMRG, 
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we use the extended infinite-size algorithm by Hida,^-* which enables us to treat non-uniform 
lattice models also. 



3. Results 

3. 1 Response of the two different insulators 

3.1.1 Mott insulator 

We calculate the polarization P in the presence of interaction for the pure systems by the 
DMRG. For V = 1.0, where the ground state of the system is the Luttinger liquid. Then the 
P — E curve is smooth as shown in Fig. l.-*^^) However for V = 6.0, where the system is in a 
Mott insulator phase , the P — E curve exhibits a stepwise behavior. Since the polarization 
corresponds to the mass center of the electrons, these steps represent a discontinuous change 
of the charge configuration. These steps actually come from the generation of kink-pairs. In 
the inset of Fig. 1, the charge distribution aX, E = E^ rii AE' are shown, where E^ is the 
smallest value of the electric field at which the kinks are generated. At the step, the kink-pair 
is generated and the electrons between the kinks are shifted by one site to the right which 
compensates for the collapse of the Mott gap. When the kinks are separated by length 
the number of electrons between kinks is 1/2. Then the energy gain of the length I kinks is 
estimated as ^ El/2. Then the first step is due to the kinks with longest length L. 

3.1.2 Anderson insulator 

We also calculated the polarization P for the Anderson insulator. Fig. 2 is the P — E 
curve for the randomness strength W = b without interaction. In the presence of randomness, 
the P — E curve exhibits a stepwise behavior also ,which is caused by crossing of one particle 
localized states. Matrix elements of the Hamiltonian between arbitrary localized states are 
exponentially small as a function of the distance. Therefore energy repulsion between the 
highest occupied and the lowest unoccupied states at both ends of the system is practically 
negligible. Then charge reconstruction of the ground state occurs which corresponds to the 
electron tunneling between the localized states. The inset of Fig. 2 is the charge distributions 
at = =t where Eg is the field strength of the step. The charge distribution near the 
center of the lattice is unchanged when E crosses the step, but that near the edges is modified. 

3.2 Linear Response Regime 

3.2.1 Susceptibility of non-interacting systems 

In this section we focus on the linear response region(£' -^t/L). We calculate the suscep- 
tibility X as a function of L by numerically differentiating the ground state energy. At first, 
we consider the non- interacting case(F = 0). In order to calculate x, we apply small electric 
field ±AE. LAE ~ lO^^. We need to avoid the occurence of the step in ±AE. When the 
step occurs near E = accidentally, we do not use the data. The number of occurence of 
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such step is typicaUy 1 2 in 100 samples when W is large(W^ ^ 5). The calculations are 



the system size is sufficiently large. Therefore we take an average of log x as is shown in Fig. 
3. We averaged over 500 realizations of the disorder potentials. For W = 0, x increases as 
X ~ and x seems to diverge. This implies the system is metallic. On the other hand, for 
W = 3, X convergent to a finite value which is consistent with the insulating ground state. 
For intermediate W, saturation of x is not clearly observed up to L = 500 and only deviation 
from X is observed. This is because the localization length ^ is larger than the system 

length we used. 

Therefore we need a finite-size scaling analysis to determine localization length ^. We 
perform the finite-size scaling analysis by assuming 



where ^{W) is a localization length and g{x) is a scaling function. If L is much larger than 
^, X becomes constant. Therefore the scaling function g behaves as g{x) ~ for x ^ 0. On 
the other hand, when L is much smaller than ^ (L <S ^), the wavefunction spread over entire 
system and the system looks like metallic. Then, % ^ is expected and the scaling function 
g behaves as g{x) ~ const, for x — > oo. The scaling function g we obtained for F = is shown 
in Fig. 3. It shows the above-mentioned one parameter scaling hypothesis works quite well in 
the model. 

3.2.2 Metal-insulator transition in attractive V 

Using the random DMRG, we study the systems with attractive electron-electron interac- 
tion. In our implementation of the DMRG, three or four finite lattice sweeps are performed to 
get the convergence of the ground state energy and the number of retained states for each block 
is 60-100 to keep the truncation error to be less than 10~^. The susceptibility x is calculated 
for various values of V and W. The calculations are carried out up to the system size L = 100. 
We take an average over 128 realizations of the disorder potentials for F = —1.4 and over 64 
realizations of the disorder potentials for V = —0.5, —0.8, —1.0, —1.1, —1.2, —1.3, —1.5, —1.6 
and —1.8. The steps in P — £^ curve become smoother when interaction is introduced. When 
log X is larger than 5, we consider this as an influence of the step and we do not use the data. 
Then we perform the finite-size scaling analysis similarly to the noninteracting case assuming 
the same one parameter scaling hypothesis. In order to avoid finite size effect, we use the 
data for L > 20 to determine the localization length ^. The obtained scaling function g for 
F = — 1.4 is shown in Fig. 4. The scaling hypothesis for the interacting case also seems to 
work well in the present model. For V = —1.8 the L dependence of x deviates from x ~ in 
L < 50 even for = 0. We could only treat systems V > —1.6 due to the finite size effects. 

In Fig. 5, W dependence of the localization length ^ is shown for various value of V. 



always carried out in the localized region except W = 0, then we expect X ~ ^ 



when 




(5) 
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Note that the localization length is normalized by ^(VF = 3). For V = 0, it seems that ^ 
is divergent at = 0, but for V = —lA, ^ increases rapidly with decreasing W and diverge 
around W ^ 1. 

In order to determine the critical disorder strength Wc, we fit the localization length ^ as 

^{W) = {A + B{W-W^)){W- W,f, (6) 

where A, B and (3 are fitting parameters. For V = 0, we obtained Wc = 0.02 zt 0.09 and 
P = —2.1 lb 0.3, which is consistent with expected value Wc = and (3 = —2. 

In Fig. 6, the obtained phase diagram is shown. The vertical fines represent errorbars. This 
phase diagram is consistent with the one obtained from the analysis of the phase sensitibity.^^^ 

4. Summary 

In the present paper we have studied the dielectric response of the one-dimensional spin- 
less fermion model with interaction and disorder by using the random DMRG. At first, we 
have calculated the polarization of the Mott insulator and the Anderson insulator. In the 
P — E curve we observed stepwise behaviors both for the Mott insulators and the Anderson 
insulators. From the change of the charge distribution, wc could understand the stepwise be- 
haviors. For the Mott insulator the steps come from the generation of kink-pair and occur 
at £^ ~ 1/L. On the other hand, the steps for the Anderson insulator represent the crossing 
of the energy levels and the steps occur ai E l/L'^. From the zero field susceptibility, we 
performed the finite-size scaling and determined the localization length ^. Also we confirmed 
the existence of the metallic region in attractive interacting regime. 

The computation in this work has been done in part using the facilities of the Supercom- 
puter Center, ISSP, University of Tokyo. Y. H. was supported in part by a Grant-in- Aid from 
the Ministry of Education, Science and Culture of Japan and also by the Kawasaki steel 21st 
Century Foundation. 
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Fig. 1. Polarization P as a function of the applied electric field E. (a) For the pure system(VF = 0). 
The Luttinger liquid regime {V = 1.0) and the Mott insulator regime {V = 6.0). The system 
size is L = 30. (b) The charge distribution at each site, solid line -.E = Eg — AE, broken line 
■.E = Es + AE where Eg is the value of the electric field where the stepwise behavior is observed. 
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Fig. 2. The polarization P as a function of the apphed electric field non-interacting system(V^ = 
0,W = 5,L = 20). Inset: the charge distribution at E = Eg + Ai?(broken line) and E = Eg — 
Ai;(solid line). 




log(1/L) log(^L) 

Fig. 3. Left: Average logx as a function of log The system sizes are between 10 and 500. V = 0. 
W = 0, 0.5, 1.5, 3.0 and 5.0. We averaged over 500 realizations of the disorder potentials. Right: 
The scaling function 5 for F = is shown. 
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Fig. 4. Scaling plots of x/-^^ vs ^/L in a log-log scale. V — —lA. 
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Fig. 5. ^ vs y = and —1.4. ^ is scaled by the value at W = 3. The dashed lines show fitting 
curves by eq (6). 
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Fig. 6. Phase diagram of the half-filled spinless fermion model. The vertical lines are errorbars. 
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